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Abstract. Electronic transport through rubrene single-crystal field effect transistors 
(FETs) is investigated experimentally in the high carrier density regime (n ~ 0.1 
carrier/molecule). In this regime, we find that the current does not increase linearly 
with the density of charge carriers, and tends to saturate. At the same time, the 
activation energy for transport unexpectedly increases with increasing n. We perform a 
theoretical analysis in terms of a well-defined microscopic model for interacting Frohlich 
polarons, that quantitatively accounts for our experimental observations. This work is 
particularly significant for our understanding of electronic transport through organic 
FETs. 
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1. Introduction 

The use of single- crystalline material for the fabrication of organic field-effect transistors 
(FETs) has given, over the past few years, the experimental control needed for 
the investigation of the intrinsic transport properties of dielectric/organic interfaces 
[Ij. This has resulted in the observation of anisotropic transport [2], a metallic-like 
temperature dependence of the mobility [3], Hall-effect [1], and quasiparticle response 
in the infrared conductivity [5J. Following this progress, the successful quantitative 
analysis of experiments in terms of a simple microscopic model has been recently possible 
in single-crystal FETs with highly polarizable gate dielectrics [6]. In these devices, 
charge carriers accumulated electrostatically at the surface of the organic semiconductor 
were shown to interact strongly with the polarization in the dielectric, leading to the 
formation of Frohlich polarons [7]. The quantitative agreement between the behavior 
predicted by the Frohlich Hamiltonian and the experimental data demonstrates that 
the microscopic processes responsible for transport in organic FETs can be qualitatively 
different than in bulk organic semiconductors. 

Here we investigate the transport properties of rubrene (C42H28) single-crystal 
FETs with highly polarizable Ta205 gate dielectrics, in the yet unexplored high carrier 
density regime. We find that in this regime the electrical characteristics of the devices 
exhibit pronounced deviations from those of conventional FETs. Specifically, in all 
devices the source-drain current Igd stops increasing linearly with the gate voltage Vg, 
and shows a clear saturation. Concomitantly, the activation energy for temperature- 
dependent transport increases when Vg increases, a trend opposite to that usually 
observed in organic transistors. We build on the Frohlich model, that was used in our 
previous study to describe the dressing of the carriers by the polarizability of the gate 
dielectric [6], and find that the observed behavior can be quantitatively explained by 
considering the effects of the Coulomb interactions between holes. Our results confirm 
that organic single- crystal transistors are suitable for the experimental investigation 
of the intrinsic electronic properties of dielectric/organic interfaces, and extend our 
fundamental understanding of transport in organic transistors to the high carrier density 
regime. 

2. Experimental results 

The rubrene transistors used in our experiments are fabricated by laminating thin (< 
Ifim thick) rubrene single-crystals onto a substrate with pre-fabricated FET circuitry, as 
described in Refs. [HI [9]. The gate dielectric is a ~ 500 nm-thick layer of Ta205 (dielectric 
constant = 25; breakdown voltage ~ Q.5MV/cm), enabling the accumulation of up 
to 10^^ holes/cm^ in the FET channel. During device operation, however, we never 
exceed gate voltage values corresponding to half the breakdown field, to minimize the 
chance of device failure [10] (we reach hole densities of 5- 10^^ cm~^, i.e. one order 
of magnitude higher than in Ref. [6]). The transistor electrical characteristics were 
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Figure 1. Source-drain current vs. gate voltage measured for two dif- 
ferent devices at different temperatures: a) from top to bottom, T = 
280, 265, 250, 235, 215, 200, 175, 150, 125 K (measurements were performed at a source- 
drain bias Vsd — —16 V; the channel width and length are W = 340/im and L = 200/xm, 
respectively); b) from top to bottom, T = 300,260,240,220,210 K {Vsd = -1 V, 
W = 28yLtm, L = 630^m). The full lines are the theoretical fits. The inset in a) 
shows the experimentally measured activation energy of Isd plotted vs. carrier density 
n = C\Vg — Vth\/&, compared with the theory (full line). The experimental error is 
comparable with the size of the data points. The device is sketched in panel b). 

measured in the vacuum (10~^ mbar) chamber of a flow cryostat, between 300 K and 
125-210 K depending on the specific device (at low temperature the devices can easily 
break due to the difference in thermal expansion between crystal and substrate). 

Figure 1 shows the source-drain current Igd measured on two different devices as 
a function of gate voltage V^, at different temperatures. Since the carrier density n in 
the channel is linearly proportional to the gate voltage Vg [n = C{Vg — Vth)/e, with 
C capacitance per unit area and Vth threshold voltage], it is normally expected that 
Isd also increases linearly with Vg. The data, however, show a pronounced deviation 
from the conventional linear behavior, since at high Vg the source-drain current tends 
to saturate. The effect is reversible and reproducible: the devices can be measured 
many times without noticeable difference in the data, which implies that the saturation 
of Isd is not due to device degradation [10]. We have measured tens of similar devices 
(room-temperature mobility values between 1.0 and 1.5 cm^/Vs) and found a similarly 
pronounced saturation of the source-drain current at high gate voltage in virtually all 
cases, albeit with some differences in the details of the Isd — Vg characteristics (illustrated 
by the data in Fig. la and b). 
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The temperature dependence of the current Isd is also unusual. At any fixed value 
of carrier density n, the current Ig^ decreases in a thermally activated way [6J. The 
inset of Fig. la shows that the activation energy Ea depends on n and that, for densities 
larger than ~ 0.02 holes/molecule, Ea increases with increasing n. This observation is 
surprising, because normally e.g. in organic thin-film FETs) Ea exhibits the opposite 
trend, which is attributed to the effect of disorder (filling of traps) in the organic 
material. 

A behavior of the Isd vs. Vg curves such as the one just described has not been 
reported earlier in FETs realized on other conventional gate dielectrics (e.g., Si02, 
AI2O3, Si3N4). Only recently, a strongly non-linear Igd — Vg relation was observed in 
organic FETs with gate electrolytes [Ullll]- Gate electrolytes enable the accumulation 
of large density of carriers, 5 x 10^^ carriers/cm^ or more, similar to our Ta205 dielectrics. 
At these densities, at which the charges accumulate in the uppermost molecular layer 
of the crystal [6] , the average distance between the carriers is only a few molecules and 
the resulting (bare) Coulomb interaction is a few hundreds of milli-electronvolts, much 
larger than the thermal energy at room temperature. This, together with the poor 
screening associated to the thermally activated motion of charge carriers, suggests that 
interactions may be responsible for the observed anomalous FET characteristics. 

3. Microscopic model 

As we proceed to show, the inclusion of Coulomb interactions between holes does 
indeed successfully account for our experimental findings. To introduce the theoretical 
framework needed for the quantitative description of the device characteristics we 
recall the conclusions of Ref. [6], namely that holes in rubrene FETs with highly 
polarizable gate dielectrics form small Frohlich polarons due to their strong interaction 
with the ionic polarization at the interface. This conclusion was established by analyzing 
transport measurements in the low density regime, performed on FETs with a range 
of different gate dielectrics. Specifically, it was shown that, using gate dielectrics of 
increasing polarizability, the "metallic-like" transport previously observed in rubrene 
[21 [13] is progressively suppressed, turning into an activated behavior characteristic of 
self-localized particles. The experimentally observed behavior was consistently explained 
in terms of the Frohlich model -that describes the polar interaction with interface 
phonon modes- by including explicitly the narrow-band nature of the organic crystal 
[6]. In FETs with Ta205 devices, the most polarizable dielectric used in our previous 
work, the interaction with the interface modes completely dominates the carrier motion, 
reducing the mobility at 200 K by two orders of magnitude with respect to isolated 
rubrene. 

Here, we extend the analysis of Ref. [6] to the high-density regime that had not 
been analyzed previously, by adding a term which accounts for the mutual Coulomb 
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interactions between the carriers [H]. The resulting Hamiltonian: 

K 1 
H = -tJ2 + c+Ci) + ^Y^^i 9ieniXi + niVijUj. (1) 

<ij> I i,e ij^j 

is the simplest model that catches the main microscopic mechanisms that determine 
the transport in our devices: finite electronic bandwidth, polar coupling of the carriers 
with the gate dielectric, and mutual interaction between carriers. The first three terms 
describe free holes in a tight-binding scheme {cf,Ci are the corresponding creation and 
annihilation operators on site i), interacting with lattice deformations through a non- 
local coupling Qii. Physically, ga represents the polar interaction between the holes in the 
two-dimensional conducting channel and the phonons at the organic/dielectric interface 
p!5| [6]. When this interaction is sufficiently strong, as is the case for the rubrene/Ta205 
interface, it leads to the formation of small polarons, i.e. self-trapped states whose 
characteristic radius is comparable with the lattice spacing. This is demonstrated 
theoretically in Appendix A, and is confirmed by the experimental data, which exhibit 
a thermally activated mobility characteristic of small polarons (see Ref. [6] and Fig. [2] 
below). As in Ref. [6], we explicitly neglect the kinetic energy of the phonons, which is 
appropriate in the adiabatic regime, for which the lattice dynamics is slow compared to 
the band electrons. 

The last term — not considered previously — represents the mutual interaction 
between holes through the unscreened Coulomb repulsion Vij, which is long-ranged 
(this contrasts with the Hubbard model, where a purely on-site interaction is 
considered, leading to the appearance of interaction effects only at densities close 
to one carrier /molecule). In the following we evaluate quantitatively the effects of 
Coulomb interactions on the transport characteristics. We shall only present the main 
formulas needed for the comparison with the experiments, leaving the full details of the 
calculations to a further publication. 

In the temperature range of the experiment the motion of small Frohlich polarons 
at the rubrene/Ta205 interface takes place through successive uncorrelated hops from 
a given molecule to a nearest neighbor. In this case, the mobility can be rigorously 
determined by analyzing the hopping rate of a given carrier between two neighboring 
molecules (say i = 1,2) [T7]. Integrating out the electronic degrees of freedom, it is found 
that the hole dynamics is coupled to the long-range lattice polarization solely through 

I It should be noted that the electron-lattice interaction in Eq. ([T|) can mediate an effective attraction 
between holes, leading under specific conditions to the stabilization of bipolarons, i.e. two holes bound 
together by a common lattice deformation. This requires in particular that the phonon-mediated 
attraction is sufficiently strong to overcome the instantaneous Coulomb repulsion between holes. This 
condition can be cast in terms of the effective dielectric constants of the interface at high and low 
frequencies, and corresponds to 77 = (eoo + f^)/{^s + f^) being less than a given critical value rjc- An 
accurate estimate 77c = 0.131 was given in Ref. |16j for a pure 2D Frohlich interaction in the continuum 
limit. Such critical value should be further reduced in the present case, because the finite distance z of 
the conducting channel to the interface suppresses the electron-phonon interaction at large momenta, 
Mq ~ e.~'^^ / ^Jq (see Appendix A). For the present rubrene/Ta205 interface we have 77 = 0.26, which 
lies safely above the critical value rj^ so that bipolaron formation can in principle be excluded. 
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the collective coordinate Q = ^Y.e{9u — g2i)Xi. Taking advantage of the adiabatic 
assumption, we end up with the following double-well potential: 

K 



+ (2) 



where we have defined g"^ = J2e 9u{.9u—g2i) ■ The physical meaning of the above equation 
is that in the adiabatic regime, the hole follows instantaneously the dynamics of the slow 
coordinate Q in the effective potential Ead- The two minima at Q — ^g/2K are then 
associated with the hole being at site 1 and 2 respectively. The corresponding hopping 
rate can be evaluated by standard techniques jl8] , by studying the escape rate of the 
coordinate Q over the barrier. The variable ^ that appears in Eq. ([2]) (defined below) 
accounts for the instantaneous repulsion of all the other carriers in the conducting layer, 
and vanishes in the low density limit. In that case the double-well potential Eq. ([2]) 
is symmetric, with a barrier given by A = g"^ /'iK — t (the polaronic gap, valid in the 
strong coupling regime g"^ /'iK ^ t), all the carriers diffuse with the same hopping rate 
Wo = (tt;s/27r)e~'^/'*'^"'", and we recover the mobility of independent polarons obtained in 
Ref. [6j, namely /ip(T) = (ea^/^T)wo (here Ug is the frequency of the interface optical 
phonons and a is the intermolecular distance). 

In the high density regime accessible with Ta205 gate dielectrics it becomes essential 
to include the Coulomb interactions that were not analyzed previously, which requires 
considering the case of a finite ^. From the model Eq. ([T]) it can be shown that: 

i = E (^2, - V^,)n,, (3) 

where Vij = 2e^/(es + i^)/Rij is the Coulomb potential, appropriately screened by the 
ionic polarization at the interface (k = 3 is the dielectric constant of rubrene). A 
finite ^ causes an energy unbalance between the initial and final hopping sites (the 
double- well potential EadiQ) becomes asymmetric) so that the hopping rate changes to 
w(^) = Woe"^/^'^'^"'". Now each carrier diffuses with a different hopping rate, determined 
by its own environment. Correspondingly, the mobility of the sample is obtained as the 
statistical average (w(^)) over all the possible values of the local electronic potential ^. 

Since the polaronic barrier sets the dominant energy scale, due to the high 
polarizability of the gate dielectric Ta205 used in our devices, we always have ,^ <^ A in 
the explored density range. As a result, correlations between subsequent hops can be 
neglected to a first approximation, which corresponds to replacing (w(^)) w((^)). 
This yields the leading term of the mobility in the presence of electron-electron 
interactions as 

^(T) = /ip(T)e-<«>/2fcsT_ 

The density dependent quantity (^) represents the average extra energy cost for hopping 
from site to site, induced by the long-range Coulomb interactions between the carriers. 
As the comparison between theory and experiments will show (see below), it is this 
dependence of n on the density (through (^)) that causes the saturation of Igd observed 
in Fig. 1. 
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In the linear response regime the calculation of (^) follows, through Eq. ([3]), from 
the statistical distribution of the occupation numbers {nj} constrainedto ui = 1 (site 1 is 
initially occupied by the hole, which then hops to site 2). Such constraint reflects the fact 
that the relaxation of the remaining carriers occurs on a time scale ~ Wq ^, much longer 
than that of the hopping event under consideration, ~ and can be neglected. In 
this case, the occupation numbers {uj} are given, by definition, by the pair distribution 
function g^'^\Rij) of the electronic system at equilibrium p^. Furthermore, due to the 
diffusive nature of their hopping motion (transport is thermally activated), the small 
polarons behave as classical interacting particles. The two-body correlations of such a 
classical fluid of charged particles are determined solely by the Coulomb interactions, 
and are therefore equivalent to those of a classical two-dimensional plasma. Observing 
that in the explored density range the mean inter-particle distance is always larger than 
the lattice periodicity, we can use for g^'^^ the known pair distribution function of the 
2D plasma in the continuous limit. Correspondingly, we can replace the discrete sum in 
Eq. ([3]) by an integral, which finally yields 

{0=nl dv [V{r + R12) - V{v)] g^^\r) = '^na'kBTF{V). (5) 

In the above equation, the universal function F(T) is an intrinsic property of the 
two-dimensional electron plasma. It depends on a single dimensionless parameter F 
that measures the ratio of the Coulomb interactions to the thermal energy. F{T) 
can be evaluated directly by using the data of g^'^^r) available from extensive Monte 
Carlo simulations [20]. In the range 1 < F < 20, we find F{r) = 1 -|- 0.85F to 
within 1% accuracy. In the present FET geometry, the plasma interaction parameter 
is given by F = 2y/Tme'^ / (kBT) / {tg + k), and can attain the value ~ 9 at the 
lowest temperatures/highest Vg, placing our devices in the range of weak to moderate 
interactions (electron crystallization is expected at F = 125). 

4. Discussion 

We are now in a position to compare the experimental data with our calculations. 
By inserting Eq.([5]) into Eq.(jl]), and using the definition of F, we obtain the explicit 
functional dependence of the mobility on density n = CiVg — Vth) /(^ and temperature. 
The Isd — Vg curves can then be calculated from Igd = — Vth\e^^^'^^'^''^'^ , with 

A = fipCVsdW/ L having the dimensions of conductance (here W and L are the width 
and the length of the channel, Vgd the applied source-drain voltage, and e the electron 
charge). To analyze the data, we fix the temperature and fit the Isd — Vg curve using 
/xp and Vth as fitting parameters. The polaron microscopic parameters are obtained by 
analyzing the temperature dependence of /ip extracted in this way. 

The values of /ip extracted for the different devices from the theory of interacting 
polarons are reported in Fig. [2] (red squares). They are close to the ones that one would 
obtain from linear fits of the Isd — Vg curves in the low density regime (the conventional 
definition of the mobility; data shown as black circles in Fig. [2]) and they exhibit the 
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Figure 2. Polaron mobility extracted from the data of Fig. la. Red squares are 
obtained by fitting the data with the theory for interacting polarons described in Sec. 
3, while black circles correspond to linear fits of the Isd — Vg curves restricted to the 
low density regime, as done in Ref . [6] . The lines are best fits to the polaronic thermally 
activated behavior, yielding respectively A — 53meV and A = 55meV (see text). The 
inset is the polaron mobility for the sample of Fig. lb (units are the same as in the 
main panel; the temperature range is too small to extract reliable values for A). 

same trend. Indeed, by comparing the temperature dependence of the fitted fip to the 
theoretical relation fip(T) = {ea'^uJs/2TikBT) exp{—A/kBT) we obtain, for the device 
of Fig. la, the following values of the polaron parameters: Us/ (27r) = 390cm~^ and 
A = 53meV. These values are comparable to those obtained by assuming a linear Igd 
vs. Vg dependence [us/ (27r) = 315cm~^ and A = 55meV]. This shows that the results of 
the interpretation based on the theory with interactions is compatible with the analysis 
performed by only looking at the low density part of the Igd — Vg curves, and that the 
effect of the interactions is determined by the carrier density only (i.e., without the need 
to include any additional parameter). For the second device (inset of Fig. [2]), a precise 
quantitative analysis of the temperature dependence of the mobility is prevented by 
the restricted temperature range of the available data, and it is not possible to extract 
reliable values for the parameters Us/i^^n) and A. Still, the fitted values of /xp and Vth 
are well defined for all the curves of Fig. lb, thus allowing for an accurate analysis of 
the density dependence of the mobility at each given temperature. 

The values of Vth obtained from the interacting theory also follow the same trend 
as those that one obtains by extrapolating the linear part of the Isd — Vg curves to zero 
current, i.e. the usual definition of Vth- The two estimates differ by a systematic offset 
of ~ 7 V, which is approximately the same for all samples and at all temperatures. Such 
a deviation does not have a particular physical significance, because the exact absolute 
value of Vth defined in the usual way does not have a precise physical meaning. 

The continuous lines in Fig. l.a and b are the results of the theoretical fits of 
the Isd — Vg curves, using the theory for interacting polarons. The theory reproduces 
the saturation of the source-drain current at high carrier density and the quantitative 
agreement with the data is remarkable, at all temperatures and in both the devices 
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Figure 3. Scaling plot of the data of Fig. l.a. (the same symbols have been used). 
The full line corresponds to the prediction of the 2D classical plasma. 



analyzed. From Eq. (jl]) we see that theory predicts an activation energy which increases 
with increasing density (via the term (^)/2). This is illustrated in the inset of Fig. 
la, where the circles are the experimental values extracted from the data and the 
continuous line is the theoretical curve computed using Eq.Q, again without any new 
free parameters. Also here the agreement is very satisfactory, and at small n we recover 
the value given by the analysis of the linear regime in terms of non-interacting polarons 



Notably, the theoretical derivation presented above implies that the density induced 
increase of the hopping activation barrier [cf. Eq. (jlj)] is a scaling function of the 
parameter F of the interacting plasma. This scaling behavior is checked in Fig. [3], where 
we use the experimental data to directly plot the quantity \og{Isd/A\Vg — Vth\)/\Vg — Vth\, 
which is proportional to the universal function -F(F) (see Eqs. (jl]) and ([5])), versus 



y\^g~ ^th\/T cx F. In terms of these variables and at sufficiently high gate voltages, the 
data do indeed tend to collapse on a single linear curve corresponding to F(T) = 1+0. 85F 
(full line in the figure). Considering the high sensitivity to details of such scaling plot, the 
agreement between theory and experiment is satisfactory It confirms the validity of the 
assumptions underlying our derivation, giving a strong indication that the anomalous 
behavior of the electrical characteristics observed at finite carrier concentrations is 
indeed due to the long-range Coulomb interactions between the carriers. 

5. Conclusion 

We have studied electronic transport through rubrene single-crystal FETs in the high- 
density regime and found new unexpected phenomena, such as a saturation of the 
current /^^ vs. Vg and an increase of the activation energy with increasing density. 
The experimental data are accurately reproduced by a theoretical analysis based on 
interacting small Frohlich polarons, which naturally extends previous studies in the 
low-density regime to include the effect of the long-range Coulomb repulsion between 
charge carriers. Remarkably, the quantitative description of the high-density regime 
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does not require the introduction of any new parameter, as the effect of the interactions 
is fully determined by the (known) density of charge carriers. Our results demonstrate 
that, by devoting sufficient effort to control the experimental systems, single- crystal 
organic field-effect transistors do allow a detailed quantitative study of the intrinsic 
transport properties of organic FETs. 
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Appendix A. Calculation of the polaron radius 

The theoretical analysis developed in the text treats the charge carriers as small Frohlich 
polarons, having dimensions comparable to the lattice spacing. To demonstrate that 
this is the case, in this Appendix we evaluate the radius of the polarons that form at a 
rubrene/Ta205 interface. To this aim we solve the adiabatic model Eq. variationally 
for a single hole in the HOMO band, and evaluate the polaron radius i?p as a function 
of the electron-phonon coupling strength. The calculations indeed show that i?p is of 
the order of the lattice spacing. 

The electron-phonon interaction Qij is obtained from the Fourier transform of the 
electron-phonon matrix element Mg = Moe~'^^/y/q [IHl E]. Here Mq = 2nhuJse'^l3/ S, 
with S the total surface of the device, and z the distance of the conducting layer to the 
polar dielectric. The parameter /3 = (eg — eoo) / (cs + ^)/ (coo + k) is a combination of the 
known dielectric constants of the gate dielectric and of the organic semiconductor, which 
determines the strength of the electron-phonon coupling. In principle, the interaction 
gij defined above must be cut off at short distances to account for the discreteness of 
the rubrene lattice. A precise prescription to carry out this procedure will be reported 
elsewhere. For the present purposes it suffices to say that the overall behavior of the 
polaron radius does not depend on the particular choice of the short-distance cutoff. 

As customary for lattice models, we introduce the dimensionless electron-phonon 
coupling 

A = f, (A.1) 

defined as the energy of a fully localized polaron in units of the half bandwidth D cct. 
The energy Ep can be written in general as: 



wher the sum spans all the lattice sites The model is solved by taking the following 
gaussian trial wavefunction: 



Ep 




(A.2) 



1 



exp(-a2i?,V2) 



(A.3) 



Current saturation and Coulomb interactions in organic single- crystal transistors 11 



with Af a normalization factor. The parameter a is obtained after minimization of the 
ground state energy. The polaron radius is then given by 

1/2 



(A.4) 



which tends to 1/a in the large polaron limit (small a, Rp ^ a). The results for Rp 
at different values of the electron-phonon coupling are summarized in Table lAll As 
expected, the polaron radius becomes comparable to the lattice spacing around A ~ 1, 
i.e. when the polaron energy is of the order of the half bandwidth. 



A 



Rp/a 



0.1 0.2 0.4 0.7 1.0 1.2 1.4 1.6 



6.0 3.7 2.3 1.6 1.2 0.9 0.5 0.4 0.3 



Table Al. Polaron radius in units of the lattice spacing a for different values of the 
dimensionless electron-phonon coupling A. 



Performing the sum in Eq. (1A.2|) with the electron-phonon coupling parameters 
appropriate to our devices {(3 = 0.1 and z = 2A), we obtain Ep = 170meV. The width 
of the HOMO band in the rubrene crystal has been evaluated through a semi-empirical 
ab initio calculation in Ref. [21] as 2D = 341meV. Such value should be considered as 
an upper bound to the actual bandwidth, which was shown in Ref. to be sizably 
reduced by the effects of molecular polarization. A further effective reduction of the 
bandwidth can arise from thermal fluctuations, as considered in Ref. [13]. From the 
above arguments and from the definition Eq. flA.l|) we conclude that A > 1 which, as 
can be seen from Table lAll corresponds to a small polaron. 
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